Aligning Deep Exposures of Sparse Fields

Note: The notebook in this repository 'Initialization.ipynb' goes over many of the basic concepts such as the setup of the environment/package installation and should be read first if you are new to HST images, DrizzlePac, or Astroquery.

Introduction

This notebook demonstrates aligning long exposures which have relatively few stars and a large number of cosmic rays. It is based on the example described in the ISR linked here (ACS ISR 2015-04: Basic Use of SExtractor Catalogs With TweakReg - I), but uses a much simpler methodology.

Rather than making use of external software (e.g. SExtractor) and going through the extra steps to create 'cosmic-ray cleaned' images for each visit, this notebook demonstrates new features in TweakReg designed to mitigate false detections.

TweakReg’s source finding task imagefind now includes parameters to exclude false detections and allow the software to more easily solve for the image offsets using matched sources lists. For example, dqbits is a list of DQ flag values to include as 'good' or to exclude as 'bad' before generating and matching source lists. For ACS/WFC, setting dqbits=-16 will mask hot pixels flagged by the instrument team, eliminating a common problem where TweakReg locks onto hot pixels and solves for the dither pattern. This can occur when users set the detection threshold value too low and hot pixels outnumber astronomical sources. Other new parameters allow selection for sharpness and roundness, which give users better control over source selection criteria and the mitigation of artifacts. More details on imagefindpars options may be found on the following webpage.

In [1]:
from astropy.table import Table
from astropy.io import fits
from astroquery.mast import Observations
from ccdproc import ImageFileCollection
import glob
import matplotlib.pyplot as plt
import os
import shutil
from drizzlepac import tweakreg
from drizzlepac import astrodrizzle
from IPython.display import Image
WARNING: AstropyDeprecationWarning: astropy.extern.six will be removed in 4.0, use the six module directly if it is still needed [astropy.extern.six]
The following task in the stsci.skypac package can be run with TEAL:
                                    skymatch                                    
The following tasks in the drizzlepac package can be run with TEAL:
    astrodrizzle       config_testbed      imagefindpars           mapreg       
       photeq            pixreplace           pixtopix            pixtosky      
  refimagefindpars       resetbits          runastrodriz          skytopix      
     tweakback            tweakreg           updatenpol

1. Download the Data

The data in this example are comprised of 4 exposures in the F814W filter, all from Visit 0X of HST program 10092. Each exposure was dithered by ~60 pixels along the y-axis in order to obtain coverage in the area of the CCD chip gap. The X and Y dithers are given in arcseconds by the POSTARG1 and POSTARG2 keywords recorded in the image header.

The following commands query MAST, download the calibrated, CTE-corrected FLC files, and place them in the same 'working' directory as this notebook.

In [2]:
# Query MAST for the F814W files.
science_list = Observations.query_criteria(proposal_id='10092', filters='F814W', obs_id='j8xi0x*')
Observations.download_products(science_list['obsid'], mrp_only=False, download_dir='./science',
                               productSubGroupDescription=['FLC'])

science_files = glob.glob('science/mastDownload/HST/*/*fits')
for im in science_files:
    root = im.split('/')[-1]
    os.rename(im, './'+root)
shutil.rmtree('science/')
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/j8xi0xs0q_flc.fits to ./science/mastDownload/HST/j8xi0xs0q/j8xi0xs0q_flc.fits ... [Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/j8xi0xs3q_flc.fits to ./science/mastDownload/HST/j8xi0xs3q/j8xi0xs3q_flc.fits ... [Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/j8xi0xs6q_flc.fits to ./science/mastDownload/HST/j8xi0xs6q/j8xi0xs6q_flc.fits ... [Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/j8xi0xsaq_flc.fits to ./science/mastDownload/HST/j8xi0xsaq/j8xi0xsaq_flc.fits ... [Done]

2. Inspect the image headers

In [3]:
collect = ImageFileCollection('./',
                              keywords=["asn_id", "detector", "filter2", "exptime", "postarg1", "postarg2"],
                              glob_include="*flc.fits", ext=0)
out_table = collect.summary
out_table
Out[3]:
Table masked=True length=4
fileasn_iddetectorfilter2exptimepostarg1postarg2
str18str9str3str5float64float64float64
j8xi0xs0q_flc.fitsJ8XI0X010WFCF814W507.00.062-3.0518
j8xi0xs3q_flc.fitsJ8XI0X010WFCF814W507.00.310.023747
j8xi0xs6q_flc.fitsJ8XI0X010WFCF814W507.0-0.0623.0518
j8xi0xsaq_flc.fitsJ8XI0X010WFCF814W507.00.1866.127347

3. TweakReg Alignment

Use TweakReg to align the FLC frames based on sources in the image. The provided input list (input_flc.list) is used to align the frames in the specified order with j8xi0xsaq_flc.fits as the reference image. The parameter conv_width specifies the convolution kernel width in pixels, with recommended values ~2x the PSF FWHM for detecting point sources in the FLC frame. For ACS/WFC & WFC3/UVIS, this parameter is typically set to 3.5 pixels and for WFC3/IR to 2.5 pixels, but the value can be increased in order to use compact objects such as small galaxies for alignment.

3a. Use 'default' parameters (Test1)

In [4]:
tweakreg.TweakReg('@input_flc.list',
                  imagefindcfg={'threshold': 100,'conv_width': 3.5},
                  shiftfile=True, outshifts='shift814_flc_test1.txt',
                  updatehdr=False, interactive=False, ylimit=0.4)
INFO:drizzlepac.util:Setting up logfile :  tweakreg.log
Setting up logfile :  tweakreg.log
INFO:drizzlepac.tweakreg:TweakReg Version 1.4.7(18-April-2018) started at: 01:11:45.346 (18/12/2019) 
TweakReg Version 1.4.7(18-April-2018) started at: 01:11:45.346 (18/12/2019) 
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.util:Version Information
Version Information
INFO:drizzlepac.util:--------------------
--------------------
INFO:drizzlepac.util:Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
INFO:drizzlepac.util:3.7.5 (default, Oct 25 2019, 10:52:18) 
3.7.5 (default, Oct 25 2019, 10:52:18) 
INFO:drizzlepac.util:numpy Version -> 1.17.4 
numpy Version -> 1.17.4 
INFO:drizzlepac.util:astropy Version -> 3.2.3 
astropy Version -> 3.2.3 
INFO:drizzlepac.util:stwcs Version -> 1.5.3 
stwcs Version -> 1.5.3 
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS common to all Processing Steps:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	exclusions :	
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakutils:j8xi0xsaq_flc.fits   
j8xi0xsaq_flc.fits   
INFO:drizzlepac.tweakutils:

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:Finding shifts for: 
Finding shifts for: 
INFO:drizzlepac.tweakreg:    j8xi0xsaq_flc.fits
    j8xi0xsaq_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs6q_flc.fits
    j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs3q_flc.fits
    j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs0q_flc.fits
    j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for finding sources for each input image:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	computesig :	True
INFO:drizzlepac.tweakreg:	conv_width :	3.5
INFO:drizzlepac.tweakreg:	dqbits :	
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	fluxmax :	None
INFO:drizzlepac.tweakreg:	fluxmin :	None
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	nsigma :	1.5
INFO:drizzlepac.tweakreg:	peakmax :	None
INFO:drizzlepac.tweakreg:	peakmin :	None
INFO:drizzlepac.tweakreg:	ratio :	1.0
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	roundhi :	1.0
INFO:drizzlepac.tweakreg:	roundlo :	-1.0
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	sharphi :	1.0
INFO:drizzlepac.tweakreg:	sharplo :	0.2
INFO:drizzlepac.tweakreg:	skysigma :	0.0
INFO:drizzlepac.tweakreg:	theta :	0.0
INFO:drizzlepac.tweakreg:	threshold :	100
INFO:drizzlepac.tweakreg:	updatehdr :	False
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	use_sharp_round :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakreg:	xyunits :	pixels
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xsaq_flc.fits':
===  Source finding for image 'j8xi0xsaq_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:11:45.524 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:11:45.524 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.082394
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:46.55 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 8449 objects.
     Found 8449 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:11:46.72 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:11:46.72 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 9.901318
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:47.697 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 8582 objects.
     Found 8582 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 17031
===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 17031
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs6q_flc.fits':
===  Source finding for image 'j8xi0xs6q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:11:48.050 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:11:48.050 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.081750
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:49.025 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 6124 objects.
     Found 6124 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:11:49.12 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:11:49.12 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.046386
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:50.054 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 6070 objects.
     Found 6070 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 12194
===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 12194
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs3q_flc.fits':
===  Source finding for image 'j8xi0xs3q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:11:50.453 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:11:50.453 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.116258
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:51.386 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 4621 objects.
     Found 4621 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:11:51.480 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:11:51.480 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.016353
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:52.359 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 4713 objects.
     Found 4713 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 9334
===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 9334
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs0q_flc.fits':
===  Source finding for image 'j8xi0xs0q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:11:52.72 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:11:52.72 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.574701
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:53.673 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 4884 objects.
     Found 4884 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:11:53.767 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:11:53.767 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.441343
INFO:drizzlepac.catalogs:###Source finding finished at: 01:11:54.660 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 5267 objects.
     Found 5267 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 10151
===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 10151
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Converting RA/Dec positions of reference sources from "j8xi0xsaq_flc.fits" to X,Y positions in reference WCS...
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:Performing alignment in the projection plane defined by the WCS
Performing alignment in the projection plane defined by the WCS
INFO:drizzlepac.tweakreg:derived from 'j8xi0xsaq_flc.fits'
derived from 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for matching sources:
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	searchrad :	1.0
INFO:drizzlepac.tweakreg:	searchunits :	arcseconds
INFO:drizzlepac.tweakreg:	see2dplot :	True
INFO:drizzlepac.tweakreg:	separation :	0.5
INFO:drizzlepac.tweakreg:	tolerance :	1.0
INFO:drizzlepac.tweakreg:	use2dhist :	True
INFO:drizzlepac.tweakreg:	xoffset :	0.0
INFO:drizzlepac.tweakreg:	yoffset :	0.0
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for fitting source lists:
INFO:drizzlepac.tweakreg:	fitgeometry :	rscale
INFO:drizzlepac.tweakreg:	labelsize :	8
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	nclip :	3
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	sigma :	3.0
INFO:drizzlepac.tweakreg:	ylimit :	0.4
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for creating a shiftfile:
INFO:drizzlepac.tweakreg:	outshifts :	shift814_flc_test1.txt
INFO:drizzlepac.tweakreg:	outwcs :	shifts_wcs.fits
INFO:drizzlepac.tweakreg:	shiftfile :	True
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs6q_flc.fits
Performing fit for: j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of -0.1504, -0.3482 with significance of 18.52 and 377 matches
Found initial X and Y shifts of -0.1504, -0.3482 with significance of 18.52 and 377 matches
INFO:drizzlepac.imgclasses:Found 106 matches for j8xi0xs6q_flc.fits...
Found 106 matches for j8xi0xs6q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: -0.1083  YSH: -0.3009    ROT: 359.9990055    SCALE: 0.999918
XSH: -0.1083  YSH: -0.3009    ROT: 359.9990055    SCALE: 0.999918
INFO:drizzlepac.imgclasses:FIT XRMS: 0.38       FIT YRMS: 0.39   
FIT XRMS: 0.38       FIT YRMS: 0.39   
INFO:drizzlepac.imgclasses:FIT RMSE: 0.54       FIT MAE: 0.46   
FIT RMSE: 0.54       FIT MAE: 0.46   
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 6.2e-06 (deg)   RMS_DEC: 4.3e-06 (deg)
RMS_RA: 6.2e-06 (deg)   RMS_DEC: 4.3e-06 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  106  objects.
Final solution based on  106  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 106
Total # points: 106
INFO:drizzlepac.tweakutils:# of points after clipping: 106
# of points after clipping: 106
INFO:drizzlepac.tweakutils:Total # points: 106
Total # points: 106
INFO:drizzlepac.tweakutils:# of points after clipping: 106
# of points after clipping: 106
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs6q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs3q_flc.fits
Performing fit for: j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of -0.4308, -0.1162 with significance of 14.84 and 292 matches
Found initial X and Y shifts of -0.4308, -0.1162 with significance of 14.84 and 292 matches
INFO:drizzlepac.imgclasses:Found 80 matches for j8xi0xs3q_flc.fits...
Found 80 matches for j8xi0xs3q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: -0.2161  YSH: -0.1794    ROT: 359.9979434    SCALE: 1.000018
XSH: -0.2161  YSH: -0.1794    ROT: 359.9979434    SCALE: 1.000018
INFO:drizzlepac.imgclasses:FIT XRMS: 0.45       FIT YRMS: 0.36   
FIT XRMS: 0.45       FIT YRMS: 0.36   
INFO:drizzlepac.imgclasses:FIT RMSE: 0.58       FIT MAE: 0.51   
FIT RMSE: 0.58       FIT MAE: 0.51   
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 6e-06 (deg)   RMS_DEC: 5.4e-06 (deg)
RMS_RA: 6e-06 (deg)   RMS_DEC: 5.4e-06 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  80  objects.
Final solution based on  80  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 80
Total # points: 80
INFO:drizzlepac.tweakutils:# of points after clipping: 80
# of points after clipping: 80
INFO:drizzlepac.tweakutils:Total # points: 80
Total # points: 80
INFO:drizzlepac.tweakutils:# of points after clipping: 80
# of points after clipping: 80
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs3q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs0q_flc.fits
Performing fit for: j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.156, -0.3262 with significance of 14.56 and 308 matches
Found initial X and Y shifts of 0.156, -0.3262 with significance of 14.56 and 308 matches
INFO:drizzlepac.imgclasses:Found 76 matches for j8xi0xs0q_flc.fits...
Found 76 matches for j8xi0xs0q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.1793  YSH: -0.3744    ROT: 359.998793    SCALE: 0.999986
XSH: 0.1793  YSH: -0.3744    ROT: 359.998793    SCALE: 0.999986
INFO:drizzlepac.imgclasses:FIT XRMS: 0.38       FIT YRMS: 0.31   
FIT XRMS: 0.38       FIT YRMS: 0.31   
INFO:drizzlepac.imgclasses:FIT RMSE: 0.49       FIT MAE: 0.38   
FIT RMSE: 0.49       FIT MAE: 0.38   
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 5.2e-06 (deg)   RMS_DEC: 4.4e-06 (deg)
RMS_RA: 5.2e-06 (deg)   RMS_DEC: 4.4e-06 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  76  objects.
Final solution based on  76  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 76
Total # points: 76
INFO:drizzlepac.tweakutils:# of points after clipping: 76
# of points after clipping: 76
INFO:drizzlepac.tweakutils:Total # points: 76
Total # points: 76
INFO:drizzlepac.tweakutils:# of points after clipping: 76
# of points after clipping: 76
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs0q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[4]
INFO:drizzlepac.tweakutils:Writing out shiftfile : shift814_flc_test1.txt
Writing out shiftfile : shift814_flc_test1.txt
INFO:drizzlepac.util:Trailer file written to:  tweakreg.log
Trailer file written to:  tweakreg.log
In [5]:
# Give the 'fit residual plots' a unique name for comparison with other tests.
residual_pngs = glob.glob("residual*png")
for png in residual_pngs: 
    path = os.path.abspath(os.path.join(os.curdir, png))
    new_path = os.path.abspath(os.path.join(os.curdir, 'test1_{}'.format(png)))
    os.rename(path, new_path)

To verify that TweakReg has found a good fit, it is important to inspect the fit residuals. Below are the dx,dy residuals for each FLC file with respect to the reference image j8xi0xsaq_flc.fits. TweakReg finds >100 matches per frame, but the RMS of the fit residuals is quite large, ~0.7 pixels, and the points are not nicely clustered around dx,dy=0.0, as expected for a good fit.

In [6]:
Image(filename='test1_residuals_j8xi0xs0q_flc.png')
Out[6]:
In [7]:
Image(filename='test1_residuals_j8xi0xs3q_flc.png')
Out[7]:
In [8]:
Image(filename='test1_residuals_j8xi0xs6q_flc.png')
Out[8]:
In [9]:
# Inspect the shift file for Test1
shift_table = Table.read('shift814_flc_test1.txt', format='ascii.no_header',
               names=['file', 'dx', 'dy', 'rot', 'scale', 'xrms', 'yrms'])
formats = ['.2f', '.2f', '.3f', '.5f', '.2f', '.2f']
for i, col in enumerate(shift_table.colnames[1: ]):
    shift_table[col].format = formats[i]
shift_table
Out[9]:
Table length=4
filedxdyrotscalexrmsyrms
str18float64float64float64float64float64float64
j8xi0xsaq_flc.fits0.000.000.0001.000000.000.00
j8xi0xs6q_flc.fits-0.11-0.30359.9990.999920.380.39
j8xi0xs3q_flc.fits-0.22-0.18359.9981.000020.450.36
j8xi0xs0q_flc.fits0.18-0.37359.9990.999990.380.31

3b. Adjust conv_width to find extended objects (Test 2)

In order for TweakReg to use compact galaxies rather than point sources for alignment, we increase the convolution kernel width parameter conv_width from 3.5 to 6.0 pixels in order to find sources with a FWHM ~3 pixels in the FLC frames.

In [10]:
tweakreg.TweakReg('@input_flc.list',
                  imagefindcfg={'threshold': 100, 'conv_width': 6.0},
                  shiftfile=True, outshifts='shift814_flc_test2.txt',
                  updatehdr=False, interactive=False, ylimit=0.4)
INFO:drizzlepac.util:Setting up logfile :  tweakreg.log
Setting up logfile :  tweakreg.log
INFO:drizzlepac.tweakreg:TweakReg Version 1.4.7(18-April-2018) started at: 01:12:04.43 (18/12/2019) 
TweakReg Version 1.4.7(18-April-2018) started at: 01:12:04.43 (18/12/2019) 
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.util:Version Information
Version Information
INFO:drizzlepac.util:--------------------
--------------------
INFO:drizzlepac.util:Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
INFO:drizzlepac.util:3.7.5 (default, Oct 25 2019, 10:52:18) 
3.7.5 (default, Oct 25 2019, 10:52:18) 
INFO:drizzlepac.util:numpy Version -> 1.17.4 
numpy Version -> 1.17.4 
INFO:drizzlepac.util:astropy Version -> 3.2.3 
astropy Version -> 3.2.3 
INFO:drizzlepac.util:stwcs Version -> 1.5.3 
stwcs Version -> 1.5.3 
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS common to all Processing Steps:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	exclusions :	
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakutils:j8xi0xsaq_flc.fits   
j8xi0xsaq_flc.fits   
INFO:drizzlepac.tweakutils:

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:Finding shifts for: 
Finding shifts for: 
INFO:drizzlepac.tweakreg:    j8xi0xsaq_flc.fits
    j8xi0xsaq_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs6q_flc.fits
    j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs3q_flc.fits
    j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs0q_flc.fits
    j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for finding sources for each input image:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	computesig :	True
INFO:drizzlepac.tweakreg:	conv_width :	6.0
INFO:drizzlepac.tweakreg:	dqbits :	
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	fluxmax :	None
INFO:drizzlepac.tweakreg:	fluxmin :	None
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	nsigma :	1.5
INFO:drizzlepac.tweakreg:	peakmax :	None
INFO:drizzlepac.tweakreg:	peakmin :	None
INFO:drizzlepac.tweakreg:	ratio :	1.0
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	roundhi :	1.0
INFO:drizzlepac.tweakreg:	roundlo :	-1.0
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	sharphi :	1.0
INFO:drizzlepac.tweakreg:	sharplo :	0.2
INFO:drizzlepac.tweakreg:	skysigma :	0.0
INFO:drizzlepac.tweakreg:	theta :	0.0
INFO:drizzlepac.tweakreg:	threshold :	100
INFO:drizzlepac.tweakreg:	updatehdr :	False
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	use_sharp_round :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakreg:	xyunits :	pixels
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xsaq_flc.fits':
===  Source finding for image 'j8xi0xsaq_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:04.598 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:04.598 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.082394
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:05.873 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 2697 objects.
     Found 2697 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:05.955 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:05.955 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 9.901318
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:07.309 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 2634 objects.
     Found 2634 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 5331
===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 5331
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs6q_flc.fits':
===  Source finding for image 'j8xi0xs6q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:07.603 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:07.603 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.081750
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:08.847 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1841 objects.
     Found 1841 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:08.925 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:08.925 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.046386
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:10.173 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1847 objects.
     Found 1847 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 3688
===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 3688
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs3q_flc.fits':
===  Source finding for image 'j8xi0xs3q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:10.406 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:10.406 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.116258
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:11.634 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1350 objects.
     Found 1350 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:11.707 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:11.707 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.016353
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:12.934 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1396 objects.
     Found 1396 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 2746
===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 2746
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs0q_flc.fits':
===  Source finding for image 'j8xi0xs0q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:13.171 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:13.171 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.574701
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:14.412 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1475 objects.
     Found 1475 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:14.48 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:14.48 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.441343
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:15.7 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1557 objects.
     Found 1557 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 3032
===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 3032
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Converting RA/Dec positions of reference sources from "j8xi0xsaq_flc.fits" to X,Y positions in reference WCS...
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:Performing alignment in the projection plane defined by the WCS
Performing alignment in the projection plane defined by the WCS
INFO:drizzlepac.tweakreg:derived from 'j8xi0xsaq_flc.fits'
derived from 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for matching sources:
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	searchrad :	1.0
INFO:drizzlepac.tweakreg:	searchunits :	arcseconds
INFO:drizzlepac.tweakreg:	see2dplot :	True
INFO:drizzlepac.tweakreg:	separation :	0.5
INFO:drizzlepac.tweakreg:	tolerance :	1.0
INFO:drizzlepac.tweakreg:	use2dhist :	True
INFO:drizzlepac.tweakreg:	xoffset :	0.0
INFO:drizzlepac.tweakreg:	yoffset :	0.0
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for fitting source lists:
INFO:drizzlepac.tweakreg:	fitgeometry :	rscale
INFO:drizzlepac.tweakreg:	labelsize :	8
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	nclip :	3
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	sigma :	3.0
INFO:drizzlepac.tweakreg:	ylimit :	0.4
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for creating a shiftfile:
INFO:drizzlepac.tweakreg:	outshifts :	shift814_flc_test2.txt
INFO:drizzlepac.tweakreg:	outwcs :	shifts_wcs.fits
INFO:drizzlepac.tweakreg:	shiftfile :	True
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs6q_flc.fits
Performing fit for: j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.1081, -0.2149 with significance of 34.01 and 86 matches
Found initial X and Y shifts of 0.1081, -0.2149 with significance of 34.01 and 86 matches
INFO:drizzlepac.imgclasses:Found 52 matches for j8xi0xs6q_flc.fits...
Found 52 matches for j8xi0xs6q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0241  YSH: -0.1148    ROT: 359.9993301    SCALE: 0.999993
XSH: 0.0241  YSH: -0.1148    ROT: 359.9993301    SCALE: 0.999993
INFO:drizzlepac.imgclasses:FIT XRMS: 0.086      FIT YRMS: 0.095  
FIT XRMS: 0.086      FIT YRMS: 0.095  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.13       FIT MAE: 0.096  
FIT RMSE: 0.13       FIT MAE: 0.096  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.5e-06 (deg)   RMS_DEC: 9.4e-07 (deg)
RMS_RA: 1.5e-06 (deg)   RMS_DEC: 9.4e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  44  objects.
Final solution based on  44  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 44
Total # points: 44
INFO:drizzlepac.tweakutils:# of points after clipping: 44
# of points after clipping: 44
INFO:drizzlepac.tweakutils:Total # points: 44
Total # points: 44
INFO:drizzlepac.tweakutils:# of points after clipping: 44
# of points after clipping: 44
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs6q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs3q_flc.fits
Performing fit for: j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of -0.02424, -0.05346 with significance of 30.59 and 71 matches
Found initial X and Y shifts of -0.02424, -0.05346 with significance of 30.59 and 71 matches
INFO:drizzlepac.imgclasses:Found 45 matches for j8xi0xs3q_flc.fits...
Found 45 matches for j8xi0xs3q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0649  YSH: -0.2451    ROT: 359.9996905    SCALE: 0.999982
XSH: 0.0649  YSH: -0.2451    ROT: 359.9996905    SCALE: 0.999982
INFO:drizzlepac.imgclasses:FIT XRMS: 0.055      FIT YRMS: 0.14   
FIT XRMS: 0.055      FIT YRMS: 0.14   
INFO:drizzlepac.imgclasses:FIT RMSE: 0.15       FIT MAE: 0.1    
FIT RMSE: 0.15       FIT MAE: 0.1    
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 2e-06 (deg)   RMS_DEC: 4.2e-07 (deg)
RMS_RA: 2e-06 (deg)   RMS_DEC: 4.2e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  37  objects.
Final solution based on  37  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs3q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs0q_flc.fits
Performing fit for: j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.03332, 0.1149 with significance of 32.52 and 72 matches
Found initial X and Y shifts of 0.03332, 0.1149 with significance of 32.52 and 72 matches
INFO:drizzlepac.imgclasses:Found 47 matches for j8xi0xs0q_flc.fits...
Found 47 matches for j8xi0xs0q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.1140  YSH: -0.3754    ROT: 6.329487741e-05    SCALE: 0.999962
XSH: 0.1140  YSH: -0.3754    ROT: 6.329487741e-05    SCALE: 0.999962
INFO:drizzlepac.imgclasses:FIT XRMS: 0.089      FIT YRMS: 0.12   
FIT XRMS: 0.089      FIT YRMS: 0.12   
INFO:drizzlepac.imgclasses:FIT RMSE: 0.15       FIT MAE: 0.11   
FIT RMSE: 0.15       FIT MAE: 0.11   
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.8e-06 (deg)   RMS_DEC: 9.4e-07 (deg)
RMS_RA: 1.8e-06 (deg)   RMS_DEC: 9.4e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  39  objects.
Final solution based on  39  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 39
Total # points: 39
INFO:drizzlepac.tweakutils:# of points after clipping: 39
# of points after clipping: 39
INFO:drizzlepac.tweakutils:Total # points: 39
Total # points: 39
INFO:drizzlepac.tweakutils:# of points after clipping: 39
# of points after clipping: 39
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs0q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[4]
INFO:drizzlepac.tweakutils:Writing out shiftfile : shift814_flc_test2.txt
Writing out shiftfile : shift814_flc_test2.txt
INFO:drizzlepac.util:Trailer file written to:  tweakreg.log
Trailer file written to:  tweakreg.log
In [11]:
# Give the 'fit residual plots' a unique name for comparison with subsequent tests.
residual_pngs = glob.glob("residual*png")
for png in residual_pngs: 
    path = os.path.abspath(os.path.join(os.curdir, png))
    new_path = os.path.abspath(os.path.join(os.curdir, 'test2_{}'.format(png)))
    os.rename(path, new_path)

TweakReg now finds matches ~50 objects per frame, and the fit for the first matched image j8xi0xs0q_flc.fits looks good, with an RMS ~0.1 pixels and with the residuals dx,dy clustered around 0.0. For the other two frames j8xi0xs3q_flc.fits and j8xi0xs6q_flc.fits, the RMS is ~0.2 pixels, and the points are not centered around dx,dy=0 pixels.

In [12]:
Image(filename='test2_residuals_j8xi0xs0q_flc.png')
Out[12]:
In [13]:
Image(filename='test2_residuals_j8xi0xs3q_flc.png')
Out[13]:
In [14]:
Image(filename='test2_residuals_j8xi0xs6q_flc.png')
Out[14]:
In [15]:
# Inspect the shift file for Test2
shift_table = Table.read('shift814_flc_test2.txt', format='ascii.no_header', 
                       names=['file', 'dx', 'dy', 'rot', 'scale', 'xrms', 'yrms'])
formats = ['.2f', '.2f', '.3f', '.5f', '.2f', '.2f']
for i, col in enumerate(shift_table.colnames[1:]):
    shift_table[col].format = formats[i]
shift_table
Out[15]:
Table length=4
filedxdyrotscalexrmsyrms
str18float64float64float64float64float64float64
j8xi0xsaq_flc.fits0.000.000.0001.000000.000.00
j8xi0xs6q_flc.fits0.02-0.11359.9990.999990.090.10
j8xi0xs3q_flc.fits0.06-0.25360.0000.999980.060.14
j8xi0xs0q_flc.fits0.11-0.380.0000.999960.090.12

3c. Exclude flagged pixels with dqbits (Test 3)

To further improve the alignment, we make use of flags in the DQ array of the FLC files. The source finding parameters in TweakReg may be modified to exclude flagged pixels when generating lists of objects in each image.

Setting the parameter dqbits=0 will consider all non-zero pixels in the DQ mask to be “bad” pixels, and the corresponding image pixels will not be used for source finding. The default value of 'None' will turn off the use of image’s DQ array for source finding. In this case, AstroDrizzle was already run by MAST on visit 0X, and cosmic-ray flags were populated in the DQ array of the FLC frames. Since the exposures within this visit were already well aligned, the 4096 flags for cosmic rays are useful for excluding false detections.

In [16]:
tweakreg.TweakReg('@input_flc.list',
                  imagefindcfg={'threshold': 100, 'conv_width': 6.0, 'dqbits': 0},
                  shiftfile=True, outshifts='shift814_flc_test3.txt',
                  updatehdr=False, interactive=False, ylimit=0.4)
INFO:drizzlepac.util:Setting up logfile :  tweakreg.log
Setting up logfile :  tweakreg.log
INFO:drizzlepac.tweakreg:TweakReg Version 1.4.7(18-April-2018) started at: 01:12:19.350 (18/12/2019) 
TweakReg Version 1.4.7(18-April-2018) started at: 01:12:19.350 (18/12/2019) 
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.util:Version Information
Version Information
INFO:drizzlepac.util:--------------------
--------------------
INFO:drizzlepac.util:Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
INFO:drizzlepac.util:3.7.5 (default, Oct 25 2019, 10:52:18) 
3.7.5 (default, Oct 25 2019, 10:52:18) 
INFO:drizzlepac.util:numpy Version -> 1.17.4 
numpy Version -> 1.17.4 
INFO:drizzlepac.util:astropy Version -> 3.2.3 
astropy Version -> 3.2.3 
INFO:drizzlepac.util:stwcs Version -> 1.5.3 
stwcs Version -> 1.5.3 
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS common to all Processing Steps:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	exclusions :	
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakutils:j8xi0xsaq_flc.fits   
j8xi0xsaq_flc.fits   
INFO:drizzlepac.tweakutils:

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:Finding shifts for: 
Finding shifts for: 
INFO:drizzlepac.tweakreg:    j8xi0xsaq_flc.fits
    j8xi0xsaq_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs6q_flc.fits
    j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs3q_flc.fits
    j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs0q_flc.fits
    j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for finding sources for each input image:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	computesig :	True
INFO:drizzlepac.tweakreg:	conv_width :	6.0
INFO:drizzlepac.tweakreg:	dqbits :	0
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	fluxmax :	None
INFO:drizzlepac.tweakreg:	fluxmin :	None
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	nsigma :	1.5
INFO:drizzlepac.tweakreg:	peakmax :	None
INFO:drizzlepac.tweakreg:	peakmin :	None
INFO:drizzlepac.tweakreg:	ratio :	1.0
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	roundhi :	1.0
INFO:drizzlepac.tweakreg:	roundlo :	-1.0
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	sharphi :	1.0
INFO:drizzlepac.tweakreg:	sharplo :	0.2
INFO:drizzlepac.tweakreg:	skysigma :	0.0
INFO:drizzlepac.tweakreg:	theta :	0.0
INFO:drizzlepac.tweakreg:	threshold :	100
INFO:drizzlepac.tweakreg:	updatehdr :	False
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	use_sharp_round :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakreg:	xyunits :	pixels
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xsaq_flc.fits':
===  Source finding for image 'j8xi0xsaq_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:19.534 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:19.534 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.082394
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:20.86 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1089 objects.
     Found 1089 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:20.956 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:20.956 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 9.901318
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:22.190 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 956 objects.
     Found 956 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 2045
===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 2045
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs6q_flc.fits':
===  Source finding for image 'j8xi0xs6q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:22.43 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:22.43 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.081750
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:23. (18/12/2019)
INFO:drizzlepac.catalogs:     Found 762 objects.
     Found 762 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:23.74 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:23.74 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.046386
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:24.978 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 707 objects.
     Found 707 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 1469
===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 1469
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs3q_flc.fits':
===  Source finding for image 'j8xi0xs3q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:25.221 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:25.221 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.116258
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:26.45 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 621 objects.
     Found 621 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:26.550 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:26.550 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.016353
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:27.776 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 541 objects.
     Found 541 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 1162
===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 1162
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs0q_flc.fits':
===  Source finding for image 'j8xi0xs0q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:27.995 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:27.995 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.574701
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:29.220 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 612 objects.
     Found 612 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:29.319 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:29.319 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.441343
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:30.54 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 608 objects.
     Found 608 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 1220
===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 1220
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Converting RA/Dec positions of reference sources from "j8xi0xsaq_flc.fits" to X,Y positions in reference WCS...
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:Performing alignment in the projection plane defined by the WCS
Performing alignment in the projection plane defined by the WCS
INFO:drizzlepac.tweakreg:derived from 'j8xi0xsaq_flc.fits'
derived from 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for matching sources:
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	searchrad :	1.0
INFO:drizzlepac.tweakreg:	searchunits :	arcseconds
INFO:drizzlepac.tweakreg:	see2dplot :	True
INFO:drizzlepac.tweakreg:	separation :	0.5
INFO:drizzlepac.tweakreg:	tolerance :	1.0
INFO:drizzlepac.tweakreg:	use2dhist :	True
INFO:drizzlepac.tweakreg:	xoffset :	0.0
INFO:drizzlepac.tweakreg:	yoffset :	0.0
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for fitting source lists:
INFO:drizzlepac.tweakreg:	fitgeometry :	rscale
INFO:drizzlepac.tweakreg:	labelsize :	8
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	nclip :	3
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	sigma :	3.0
INFO:drizzlepac.tweakreg:	ylimit :	0.4
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for creating a shiftfile:
INFO:drizzlepac.tweakreg:	outshifts :	shift814_flc_test3.txt
INFO:drizzlepac.tweakreg:	outwcs :	shifts_wcs.fits
INFO:drizzlepac.tweakreg:	shiftfile :	True
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs6q_flc.fits
Performing fit for: j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.2628, -0.1438 with significance of 41.42 and 56 matches
Found initial X and Y shifts of 0.2628, -0.1438 with significance of 41.42 and 56 matches
INFO:drizzlepac.imgclasses:Found 49 matches for j8xi0xs6q_flc.fits...
Found 49 matches for j8xi0xs6q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0169  YSH: -0.0960    ROT: 359.9991133    SCALE: 0.999992
XSH: 0.0169  YSH: -0.0960    ROT: 359.9991133    SCALE: 0.999992
INFO:drizzlepac.imgclasses:FIT XRMS: 0.086      FIT YRMS: 0.072  
FIT XRMS: 0.086      FIT YRMS: 0.072  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.11       FIT MAE: 0.087  
FIT RMSE: 0.11       FIT MAE: 0.087  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.2e-06 (deg)   RMS_DEC: 1e-06 (deg)
RMS_RA: 1.2e-06 (deg)   RMS_DEC: 1e-06 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  42  objects.
Final solution based on  42  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 42
Total # points: 42
INFO:drizzlepac.tweakutils:# of points after clipping: 42
# of points after clipping: 42
INFO:drizzlepac.tweakutils:Total # points: 42
Total # points: 42
INFO:drizzlepac.tweakutils:# of points after clipping: 42
# of points after clipping: 42
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs6q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs3q_flc.fits
Performing fit for: j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.02802, -0.004584 with significance of 36.42 and 51 matches
Found initial X and Y shifts of 0.02802, -0.004584 with significance of 36.42 and 51 matches
INFO:drizzlepac.imgclasses:Found 43 matches for j8xi0xs3q_flc.fits...
Found 43 matches for j8xi0xs3q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0565  YSH: -0.2828    ROT: 359.9987159    SCALE: 0.999975
XSH: 0.0565  YSH: -0.2828    ROT: 359.9987159    SCALE: 0.999975
INFO:drizzlepac.imgclasses:FIT XRMS: 0.051      FIT YRMS: 0.074  
FIT XRMS: 0.051      FIT YRMS: 0.074  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.09       FIT MAE: 0.072  
FIT RMSE: 0.09       FIT MAE: 0.072  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.1e-06 (deg)   RMS_DEC: 5.1e-07 (deg)
RMS_RA: 1.1e-06 (deg)   RMS_DEC: 5.1e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  36  objects.
Final solution based on  36  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 36
Total # points: 36
INFO:drizzlepac.tweakutils:# of points after clipping: 36
# of points after clipping: 36
INFO:drizzlepac.tweakutils:Total # points: 36
Total # points: 36
INFO:drizzlepac.tweakutils:# of points after clipping: 36
# of points after clipping: 36
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs3q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[4]
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs0q_flc.fits
Performing fit for: j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0, 0 with significance of 37.16 and 50 matches
Found initial X and Y shifts of 0, 0 with significance of 37.16 and 50 matches
INFO:drizzlepac.imgclasses:Found 44 matches for j8xi0xs0q_flc.fits...
Found 44 matches for j8xi0xs0q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.1048  YSH: -0.3954    ROT: 359.9994651    SCALE: 0.999971
XSH: 0.1048  YSH: -0.3954    ROT: 359.9994651    SCALE: 0.999971
INFO:drizzlepac.imgclasses:FIT XRMS: 0.076      FIT YRMS: 0.082  
FIT XRMS: 0.076      FIT YRMS: 0.082  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.11       FIT MAE: 0.081  
FIT RMSE: 0.11       FIT MAE: 0.081  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.3e-06 (deg)   RMS_DEC: 8.5e-07 (deg)
RMS_RA: 1.3e-06 (deg)   RMS_DEC: 8.5e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  37  objects.
Final solution based on  37  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs0q_flc.fits...
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[1]
INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[4]
INFO:drizzlepac.tweakutils:Writing out shiftfile : shift814_flc_test3.txt
Writing out shiftfile : shift814_flc_test3.txt
INFO:drizzlepac.util:Trailer file written to:  tweakreg.log
Trailer file written to:  tweakreg.log
In [17]:
# Give the 'fit residual plots' a unique name for comparison with other tests.
residual_pngs = glob.glob("residual*png")
for png in residual_pngs: 
    path = os.path.abspath(os.path.join(os.curdir, png))
    new_path = os.path.abspath(os.path.join(os.curdir, 'test3_{}'.format(png)))
    os.rename(path, new_path)   

In this third test, TweakReg again finds ~40 matches per frame, but with spurious detections eliminated has an easier time locking onto the correct solution. The fit residuals below all have an RMS ~0.1 pixels and the points are all clustered around dx,dy=0.

In [18]:
Image(filename='test3_residuals_j8xi0xs0q_flc.png')
Out[18]:
In [19]:
Image(filename='test3_residuals_j8xi0xs3q_flc.png')
Out[19]:
In [20]:
Image(filename='test3_residuals_j8xi0xs6q_flc.png')
Out[20]:
In [21]:
# Inspect the shift file for Test2
shift_table = Table.read('shift814_flc_test3.txt', format='ascii.no_header',
                         names=['file', 'dx', 'dy', 'rot', 'scale', 'xrms', 'yrms'])
formats = ['.2f', '.2f', '.3f', '.5f', '.2f', '.2f']
for i, col in enumerate(shift_table.colnames[1:]):
    shift_table[col].format = formats[i]
shift_table
Out[21]:
Table length=4
filedxdyrotscalexrmsyrms
str18float64float64float64float64float64float64
j8xi0xsaq_flc.fits0.000.000.0001.000000.000.00
j8xi0xs6q_flc.fits0.02-0.10359.9990.999990.090.07
j8xi0xs3q_flc.fits0.06-0.28359.9990.999970.050.07
j8xi0xs0q_flc.fits0.10-0.40359.9990.999970.080.08

3d. Rerun TweakReg to update the header WCS

Now run TweakReg with updatehdr=True to update the image headers with this solution.

In [22]:
tweakreg.TweakReg('@input_flc.list',
                  imagefindcfg={'threshold': 100, 'conv_width': 6.0, 'dqbits': 0},
                  shiftfile=False, updatehdr=True, interactive=False)
INFO:drizzlepac.util:Setting up logfile :  tweakreg.log
Setting up logfile :  tweakreg.log
INFO:drizzlepac.tweakreg:TweakReg Version 1.4.7(18-April-2018) started at: 01:12:33.767 (18/12/2019) 
TweakReg Version 1.4.7(18-April-2018) started at: 01:12:33.767 (18/12/2019) 
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.util:Version Information
Version Information
INFO:drizzlepac.util:--------------------
--------------------
INFO:drizzlepac.util:Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
INFO:drizzlepac.util:3.7.5 (default, Oct 25 2019, 10:52:18) 
3.7.5 (default, Oct 25 2019, 10:52:18) 
INFO:drizzlepac.util:numpy Version -> 1.17.4 
numpy Version -> 1.17.4 
INFO:drizzlepac.util:astropy Version -> 3.2.3 
astropy Version -> 3.2.3 
INFO:drizzlepac.util:stwcs Version -> 1.5.3 
stwcs Version -> 1.5.3 
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS common to all Processing Steps:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	exclusions :	
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakutils:j8xi0xsaq_flc.fits   
j8xi0xsaq_flc.fits   
INFO:drizzlepac.tweakutils:

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:Finding shifts for: 
Finding shifts for: 
INFO:drizzlepac.tweakreg:    j8xi0xsaq_flc.fits
    j8xi0xsaq_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs6q_flc.fits
    j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs3q_flc.fits
    j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:    j8xi0xs0q_flc.fits
    j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for finding sources for each input image:
INFO:drizzlepac.tweakreg:	clean :	False
INFO:drizzlepac.tweakreg:	computesig :	True
INFO:drizzlepac.tweakreg:	conv_width :	6.0
INFO:drizzlepac.tweakreg:	dqbits :	0
INFO:drizzlepac.tweakreg:	enforce_user_order :	True
INFO:drizzlepac.tweakreg:	expand_refcat :	False
INFO:drizzlepac.tweakreg:	fluxmax :	None
INFO:drizzlepac.tweakreg:	fluxmin :	None
INFO:drizzlepac.tweakreg:	input :	@input_flc.list
INFO:drizzlepac.tweakreg:	interactive :	False
INFO:drizzlepac.tweakreg:	nsigma :	1.5
INFO:drizzlepac.tweakreg:	peakmax :	None
INFO:drizzlepac.tweakreg:	peakmin :	None
INFO:drizzlepac.tweakreg:	ratio :	1.0
INFO:drizzlepac.tweakreg:	refimage :	
INFO:drizzlepac.tweakreg:	roundhi :	1.0
INFO:drizzlepac.tweakreg:	roundlo :	-1.0
INFO:drizzlepac.tweakreg:	runfile :	tweakreg.log
INFO:drizzlepac.tweakreg:	sharphi :	1.0
INFO:drizzlepac.tweakreg:	sharplo :	0.2
INFO:drizzlepac.tweakreg:	skysigma :	0.0
INFO:drizzlepac.tweakreg:	theta :	0.0
INFO:drizzlepac.tweakreg:	threshold :	100
INFO:drizzlepac.tweakreg:	updatehdr :	True
INFO:drizzlepac.tweakreg:	updatewcs :	False
INFO:drizzlepac.tweakreg:	use_sharp_round :	False
INFO:drizzlepac.tweakreg:	verbose :	False
INFO:drizzlepac.tweakreg:	writecat :	True
INFO:drizzlepac.tweakreg:	xyunits :	pixels
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xsaq_flc.fits':
===  Source finding for image 'j8xi0xsaq_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:34.128 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 1) started at: 01:12:34.128 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.082394
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:35.368 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 1089 objects.
     Found 1089 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:35.452 (18/12/2019)
  #  Source finding for 'j8xi0xsaq_flc.fits', EXT=('SCI', 2) started at: 01:12:35.452 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 9.901318
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:36.781 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 956 objects.
     Found 956 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 2045
===  FINAL number of objects in image 'j8xi0xsaq_flc.fits': 2045
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs6q_flc.fits':
===  Source finding for image 'j8xi0xs6q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:37.115 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 1) started at: 01:12:37.115 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.081750
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:38.347 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 762 objects.
     Found 762 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:38.430 (18/12/2019)
  #  Source finding for 'j8xi0xs6q_flc.fits', EXT=('SCI', 2) started at: 01:12:38.430 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.046386
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:39.655 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 707 objects.
     Found 707 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 1469
===  FINAL number of objects in image 'j8xi0xs6q_flc.fits': 1469
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs3q_flc.fits':
===  Source finding for image 'j8xi0xs3q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:39.974 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 1) started at: 01:12:39.974 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.116258
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:41.206 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 621 objects.
     Found 621 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:41.289 (18/12/2019)
  #  Source finding for 'j8xi0xs3q_flc.fits', EXT=('SCI', 2) started at: 01:12:41.289 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.016353
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:42.516 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 541 objects.
     Found 541 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 1162
===  FINAL number of objects in image 'j8xi0xs3q_flc.fits': 1162
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:===  Source finding for image 'j8xi0xs0q_flc.fits':
===  Source finding for image 'j8xi0xs0q_flc.fits':
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:42.828 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 1) started at: 01:12:42.828 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.574701
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:44.069 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 612 objects.
     Found 612 objects.
INFO:drizzlepac.catalogs:  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:44.162 (18/12/2019)
  #  Source finding for 'j8xi0xs0q_flc.fits', EXT=('SCI', 2) started at: 01:12:44.162 (18/12/2019)
INFO:drizzlepac.catalogs:   Finding sources using sky sigma = 10.441343
INFO:drizzlepac.catalogs:###Source finding finished at: 01:12:45.416 (18/12/2019)
INFO:drizzlepac.catalogs:     Found 608 objects.
     Found 608 objects.
INFO:drizzlepac.imgclasses:===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 1220
===  FINAL number of objects in image 'j8xi0xs0q_flc.fits': 1220
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Converting RA/Dec positions of reference sources from "j8xi0xsaq_flc.fits" to X,Y positions in reference WCS...
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:Performing alignment in the projection plane defined by the WCS
Performing alignment in the projection plane defined by the WCS
INFO:drizzlepac.tweakreg:derived from 'j8xi0xsaq_flc.fits'
derived from 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.tweakreg:===============================================================
===============================================================
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for matching sources:
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	searchrad :	1.0
INFO:drizzlepac.tweakreg:	searchunits :	arcseconds
INFO:drizzlepac.tweakreg:	see2dplot :	True
INFO:drizzlepac.tweakreg:	separation :	0.5
INFO:drizzlepac.tweakreg:	tolerance :	1.0
INFO:drizzlepac.tweakreg:	use2dhist :	True
INFO:drizzlepac.tweakreg:	xoffset :	0.0
INFO:drizzlepac.tweakreg:	yoffset :	0.0
INFO:drizzlepac.tweakreg:
INFO:drizzlepac.tweakreg:USER INPUT PARAMETERS for fitting source lists:
INFO:drizzlepac.tweakreg:	fitgeometry :	rscale
INFO:drizzlepac.tweakreg:	labelsize :	8
INFO:drizzlepac.tweakreg:	minobj :	15
INFO:drizzlepac.tweakreg:	nclip :	3
INFO:drizzlepac.tweakreg:	residplot :	both
INFO:drizzlepac.tweakreg:	sigma :	3.0
INFO:drizzlepac.tweakreg:	ylimit :	None
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs6q_flc.fits
Performing fit for: j8xi0xs6q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs6q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.2628, -0.1438 with significance of 41.42 and 56 matches
Found initial X and Y shifts of 0.2628, -0.1438 with significance of 41.42 and 56 matches
INFO:drizzlepac.imgclasses:Found 49 matches for j8xi0xs6q_flc.fits...
Found 49 matches for j8xi0xs6q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
Computed  rscale  fit for  j8xi0xs6q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0169  YSH: -0.0960    ROT: 359.9991133    SCALE: 0.999992
XSH: 0.0169  YSH: -0.0960    ROT: 359.9991133    SCALE: 0.999992
INFO:drizzlepac.imgclasses:FIT XRMS: 0.086      FIT YRMS: 0.072  
FIT XRMS: 0.086      FIT YRMS: 0.072  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.11       FIT MAE: 0.087  
FIT RMSE: 0.11       FIT MAE: 0.087  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.2e-06 (deg)   RMS_DEC: 1e-06 (deg)
RMS_RA: 1.2e-06 (deg)   RMS_DEC: 1e-06 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  42  objects.
Final solution based on  42  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs6q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 42
Total # points: 42
INFO:drizzlepac.tweakutils:# of points after clipping: 42
# of points after clipping: 42
INFO:drizzlepac.tweakutils:Total # points: 42
Total # points: 42
INFO:drizzlepac.tweakutils:# of points after clipping: 42
# of points after clipping: 42
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs6q_flc.fits...
....Updating header for j8xi0xs6q_flc.fits...
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',1]
Processing j8xi0xs6q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[1]
Updating header for j8xi0xs6q_flc.fits[1]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.80350213592781e-06 1.3859833511676653e-05
CD_11  CD_12: 2.80350213592781e-06 1.3859833511676653e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.3557424320416147e-05 -1.8921997422942406e-06
CD_21  CD_22: 1.3557424320416147e-05 -1.8921997422942406e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.18741245943534 2.6787674492582214
CRVAL    : 150.18741245943534 2.6787674492582214
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.04983893388069146
Plate Scale : 0.04983893388069146
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.77418556707268
ORIENTAT : 97.77418556707268
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs6q_flc.fits['SCI',2]
Processing j8xi0xs6q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs6q_flc.fits[4]
Updating header for j8xi0xs6q_flc.fits[4]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.9501899805599017e-06 1.3411066038082081e-05
CD_11  CD_12: 2.9501899805599017e-06 1.3411066038082081e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.3360383982723292e-05 -1.7270171766689415e-06
CD_21  CD_22: 1.3360383982723292e-05 -1.7270171766689415e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.21603635081973 2.674942533790987
CRVAL    : 150.21603635081973 2.674942533790987
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.049255659028792816
Plate Scale : 0.049255659028792816
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.33791051267566
ORIENTAT : 97.33791051267566
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.imgclasses:    Writing out new WCS to alternate WCS: "C"
INFO:drizzlepac.imgclasses:Updating WCSCORR table with new WCS solution "TWEAK"
WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs3q_flc.fits
Performing fit for: j8xi0xs3q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs3q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0.02802, -0.004584 with significance of 36.42 and 51 matches
Found initial X and Y shifts of 0.02802, -0.004584 with significance of 36.42 and 51 matches
INFO:drizzlepac.imgclasses:Found 43 matches for j8xi0xs3q_flc.fits...
Found 43 matches for j8xi0xs3q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
Computed  rscale  fit for  j8xi0xs3q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.0565  YSH: -0.2828    ROT: 359.9987159    SCALE: 0.999975
XSH: 0.0565  YSH: -0.2828    ROT: 359.9987159    SCALE: 0.999975
INFO:drizzlepac.imgclasses:FIT XRMS: 0.051      FIT YRMS: 0.074  
FIT XRMS: 0.051      FIT YRMS: 0.074  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.09       FIT MAE: 0.072  
FIT RMSE: 0.09       FIT MAE: 0.072  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.1e-06 (deg)   RMS_DEC: 5.1e-07 (deg)
RMS_RA: 1.1e-06 (deg)   RMS_DEC: 5.1e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  36  objects.
Final solution based on  36  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs3q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 36
Total # points: 36
INFO:drizzlepac.tweakutils:# of points after clipping: 36
# of points after clipping: 36
INFO:drizzlepac.tweakutils:Total # points: 36
Total # points: 36
INFO:drizzlepac.tweakutils:# of points after clipping: 36
# of points after clipping: 36
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs3q_flc.fits...
....Updating header for j8xi0xs3q_flc.fits...
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',1]
Processing j8xi0xs3q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[1]
Updating header for j8xi0xs3q_flc.fits[1]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.8035571769654962e-06 1.386057496912875e-05
CD_11  CD_12: 2.8035571769654962e-06 1.386057496912875e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.355770237160731e-05 -1.8925177157864187e-06
CD_21  CD_22: 1.355770237160731e-05 -1.8925177157864187e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.18823614067068 2.6785547459254855
CRVAL    : 150.18823614067068 2.6785547459254855
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.04983908423250734
Plate Scale : 0.04983908423250734
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.7750651420869
ORIENTAT : 97.7750651420869
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs3q_flc.fits['SCI',2]
Processing j8xi0xs3q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs3q_flc.fits[4]
Updating header for j8xi0xs3q_flc.fits[4]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.9502413433048623e-06 1.3411773579244721e-05
CD_11  CD_12: 2.9502413433048623e-06 1.3411773579244721e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.336066410079676e-05 -1.727323837892086e-06
CD_21  CD_22: 1.336066410079676e-05 -1.727323837892086e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.21686046607624 2.6747297830554992
CRVAL    : 150.21686046607624 2.6747297830554992
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.04925582375646241
Plate Scale : 0.04925582375646241
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.33881631845084
ORIENTAT : 97.33881631845084
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.imgclasses:    Writing out new WCS to alternate WCS: "C"
INFO:drizzlepac.imgclasses:Updating WCSCORR table with new WCS solution "TWEAK"
WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

INFO:drizzlepac.tweakreg:

INFO:drizzlepac.tweakreg:====================
====================
INFO:drizzlepac.tweakreg:Performing fit for: j8xi0xs0q_flc.fits
Performing fit for: j8xi0xs0q_flc.fits
INFO:drizzlepac.tweakreg:

INFO:drizzlepac.imgclasses:Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
Matching sources from 'j8xi0xs0q_flc.fits' with sources from reference image 'j8xi0xsaq_flc.fits'
INFO:drizzlepac.imgclasses:Computing initial guess for X and Y shifts...
Computing initial guess for X and Y shifts...
INFO:drizzlepac.imgclasses:Found initial X and Y shifts of 0, 0 with significance of 37.16 and 50 matches
Found initial X and Y shifts of 0, 0 with significance of 37.16 and 50 matches
INFO:drizzlepac.imgclasses:Found 44 matches for j8xi0xs0q_flc.fits...
Found 44 matches for j8xi0xs0q_flc.fits...
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.linearfit:Performing "rscale" fit
INFO:drizzlepac.imgclasses:Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
Computed  rscale  fit for  j8xi0xs0q_flc.fits : 
INFO:drizzlepac.imgclasses:XSH: 0.1048  YSH: -0.3954    ROT: 359.9994651    SCALE: 0.999971
XSH: 0.1048  YSH: -0.3954    ROT: 359.9994651    SCALE: 0.999971
INFO:drizzlepac.imgclasses:FIT XRMS: 0.076      FIT YRMS: 0.082  
FIT XRMS: 0.076      FIT YRMS: 0.082  
INFO:drizzlepac.imgclasses:FIT RMSE: 0.11       FIT MAE: 0.081  
FIT RMSE: 0.11       FIT MAE: 0.081  
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:RMS_RA: 1.3e-06 (deg)   RMS_DEC: 8.5e-07 (deg)
RMS_RA: 1.3e-06 (deg)   RMS_DEC: 8.5e-07 (deg)
INFO:drizzlepac.imgclasses:

INFO:drizzlepac.imgclasses:Final solution based on  37  objects.
Final solution based on  37  objects.
INFO:drizzlepac.imgclasses:Creating catalog for the fit: j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
wrote XY data to:  j8xi0xs0q_flc_catalog_fit.match
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.tweakutils:Total # points: 37
Total # points: 37
INFO:drizzlepac.tweakutils:# of points after clipping: 37
# of points after clipping: 37
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:....Updating header for j8xi0xs0q_flc.fits...
....Updating header for j8xi0xs0q_flc.fits...
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',1]
Processing j8xi0xs0q_flc.fits['SCI',1]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[1]
Updating header for j8xi0xs0q_flc.fits[1]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.803777543564549e-06 1.3861066997825172e-05
CD_11  CD_12: 2.803777543564549e-06 1.3861066997825172e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.3558173480915093e-05 -1.892639944480594e-06
CD_21  CD_22: 1.3558173480915093e-05 -1.892639944480594e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.18909423956995 2.678510113120627
CRVAL    : 150.18909423956995 2.678510113120627
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.049840714278485104
Plate Scale : 0.049840714278485104
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.77528851966095
ORIENTAT : 97.77528851966095
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Processing j8xi0xs0q_flc.fits['SCI',2]
Processing j8xi0xs0q_flc.fits['SCI',2]
INFO:drizzlepac.updatehdr:

INFO:drizzlepac.updatehdr:Updating header for j8xi0xs0q_flc.fits[4]
Updating header for j8xi0xs0q_flc.fits[4]
INFO:drizzlepac.updatehdr:    with WCS of
INFO:stwcs.wcsutil.hstwcs:WCS Keywords
WCS Keywords
INFO:stwcs.wcsutil.hstwcs:

INFO:stwcs.wcsutil.hstwcs:CD_11  CD_12: 2.9504689632131873e-06 1.3412263079601292e-05
CD_11  CD_12: 2.9504689632131873e-06 1.3412263079601292e-05
INFO:stwcs.wcsutil.hstwcs:CD_21  CD_22: 1.3361125620009374e-05 -1.7274336559109079e-06
CD_21  CD_22: 1.3361125620009374e-05 -1.7274336559109079e-06
INFO:stwcs.wcsutil.hstwcs:CRVAL    : 150.2177186165875 2.674684887108254
CRVAL    : 150.2177186165875 2.674684887108254
INFO:stwcs.wcsutil.hstwcs:CRPIX    : 2048.0 1024.0
CRPIX    : 2048.0 1024.0
INFO:stwcs.wcsutil.hstwcs:NAXIS    : 4096 2048
NAXIS    : 4096 2048
INFO:stwcs.wcsutil.hstwcs:Plate Scale : 0.049257433630216355
Plate Scale : 0.049257433630216355
INFO:stwcs.wcsutil.hstwcs:ORIENTAT : 97.33901287372346
ORIENTAT : 97.33901287372346
INFO:drizzlepac.updatehdr:WCSNAME  :  TWEAK
WCSNAME  :  TWEAK
INFO:drizzlepac.imgclasses:    Writing out new WCS to alternate WCS: "C"
INFO:drizzlepac.imgclasses:Updating WCSCORR table with new WCS solution "TWEAK"
WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

INFO:drizzlepac.imgclasses:    Saving Primary WCS to alternate WCS: "B"
INFO:drizzlepac.imgclasses:Updating WCSCORR table with new WCS solution "TWEAK"
WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

WARNING:py.warnings:/Users/etollerud/miniconda3/envs/nb-build/lib/python3.7/site-packages/astropy/io/fits/card.py:267: VerifyWarning: Keyword name 'IDCSCALEC' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created.
  keyword), VerifyWarning)

INFO:drizzlepac.util:Trailer file written to:  tweakreg.log
Trailer file written to:  tweakreg.log

4. Drizzle the aligned frames

Combine the aligned FLC files with AstroDrizzle. The ACS team now corrects for stable hot pixels (DQ flag=16) via the dark reference files, so these pixels can be considered 'good'. Full well saturated pixels (DQ flag=256) and warm pixels (DQ flag=64) may also be treated as good. More details on the recommended drizzle parameters for ACS may be found in ISR 2017-02.

In [23]:
astrodrizzle.AstroDrizzle('*flc.fits',
                          output='f814w',
                          preserve=False,
                          clean=True,
                          build=False,
                          context=False,
                          driz_sep_bits='256,64,16',
                          combine_type='median',
                          final_bits='256,64,16',
                          runfile='f814w_driz.log')
INFO:drizzlepac.util:Setting up logfile :  f814w_driz.log
Setting up logfile :  f814w_driz.log
INFO:drizzlepac.astrodrizzle:AstroDrizzle Version 3.1.3 (2019-11-06 14:37:50 -0500) started at: 01:12:51.881 (18/12/2019)
AstroDrizzle Version 3.1.3 (2019-11-06 14:37:50 -0500) started at: 01:12:51.881 (18/12/2019)
INFO:drizzlepac.astrodrizzle:

INFO:drizzlepac.astrodrizzle:Version Information
INFO:drizzlepac.astrodrizzle:--------------------
INFO:drizzlepac.astrodrizzle:Python Version [Clang 4.0.1 (tags/RELEASE_401/final)]
INFO:drizzlepac.astrodrizzle:3.7.5 (default, Oct 25 2019, 10:52:18) 
INFO:drizzlepac.astrodrizzle:numpy Version -> 1.17.4 
INFO:drizzlepac.astrodrizzle:astropy Version -> 3.2.3 
INFO:drizzlepac.astrodrizzle:stwcs Version -> 1.5.3 
INFO:drizzlepac.util:==== Processing Step  Initialization  started at  01:12:51.88 (18/12/2019)
==== Processing Step  Initialization  started at  01:12:51.88 (18/12/2019)
INFO:drizzlepac.processInput:Executing serially
INFO:drizzlepac.processInput:Setting up output name: f814w_drc.fits
INFO:drizzlepac.processInput:-Creating imageObject List as input for processing steps.
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 55.52809143066406
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 55.52809143066406
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.70233535766601
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.70233535766601
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.87463760375977
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.87463760375977
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.27589797973633
INFO:drizzlepac.imageObject:Reading in MDRIZSKY of 50.27589797973633
INFO:drizzlepac.imageObject:No context image will be created for j8xi0xs0q_flc.fits
INFO:drizzlepac.imageObject:No context image will be created for j8xi0xs3q_flc.fits
INFO:drizzlepac.imageObject:No context image will be created for j8xi0xs6q_flc.fits
INFO:drizzlepac.imageObject:No context image will be created for j8xi0xsaq_flc.fits
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs0q_flc.fits[DQ,1]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs0q_flc.fits[DQ,2]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs3q_flc.fits[DQ,1]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs3q_flc.fits[DQ,2]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs6q_flc.fits[DQ,1]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xs6q_flc.fits[DQ,2]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xsaq_flc.fits[DQ,1]
INFO:drizzlepac.resetbits:Reset bit values of 4096 to a value of 0 in j8xi0xsaq_flc.fits[DQ,2]
INFO:drizzlepac.processInput:-Creating output WCS.
INFO:drizzlepac.imageObject:No context image will be created for f814w_drc.fits
INFO:astropy.wcs.wcs:WCS Keywords
WCS Keywords
INFO:astropy.wcs.wcs:

INFO:astropy.wcs.wcs:Number of WCS axes: 2
Number of WCS axes: 2
INFO:astropy.wcs.wcs:CTYPE : 'RA---TAN'  'DEC--TAN'  
CTYPE : 'RA---TAN'  'DEC--TAN'  
INFO:astropy.wcs.wcs:CRVAL : 150.20204715246223  2.6769864498720612  
CRVAL : 150.20204715246223  2.6769864498720612  
INFO:astropy.wcs.wcs:CRPIX : 2111.5  2210.0  
CRPIX : 2111.5  2210.0  
INFO:astropy.wcs.wcs:CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
INFO:astropy.wcs.wcs:CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
INFO:astropy.wcs.wcs:NAXIS : 4223  4420
NAXIS : 4223  4420
INFO:drizzlepac.processInput:********************************************************************************
********************************************************************************
INFO:drizzlepac.processInput:*
*
INFO:drizzlepac.processInput:*  Estimated memory usage:  up to 2220 Mb.
*  Estimated memory usage:  up to 2220 Mb.
INFO:drizzlepac.processInput:*  Output image size:       4223 X 4420 pixels. 
*  Output image size:       4223 X 4420 pixels. 
INFO:drizzlepac.processInput:*  Output image file:       ~ 213 Mb. 
*  Output image file:       ~ 213 Mb. 
INFO:drizzlepac.processInput:*  Cores available:         8
*  Cores available:         8
INFO:drizzlepac.processInput:*
*
INFO:drizzlepac.processInput:********************************************************************************
********************************************************************************
INFO:drizzlepac.util:==== Processing Step  Initialization  finished at  01:12:53.699 (18/12/2019)
==== Processing Step  Initialization  finished at  01:12:53.699 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.astrodrizzle:USER INPUT PARAMETERS common to all Processing Steps:
INFO:drizzlepac.astrodrizzle:	build :	False
INFO:drizzlepac.astrodrizzle:	coeffs :	True
INFO:drizzlepac.astrodrizzle:	context :	False
INFO:drizzlepac.astrodrizzle:	crbit :	4096
INFO:drizzlepac.astrodrizzle:	group :	
INFO:drizzlepac.astrodrizzle:	in_memory :	False
INFO:drizzlepac.astrodrizzle:	input :	*flc.fits
INFO:drizzlepac.astrodrizzle:	mdriztab :	False
INFO:drizzlepac.astrodrizzle:	num_cores :	None
INFO:drizzlepac.astrodrizzle:	output :	f814w
INFO:drizzlepac.astrodrizzle:	proc_unit :	native
INFO:drizzlepac.astrodrizzle:	resetbits :	4096
INFO:drizzlepac.astrodrizzle:	runfile :	f814w_driz.log
INFO:drizzlepac.astrodrizzle:	stepsize :	10
INFO:drizzlepac.astrodrizzle:	updatewcs :	False
INFO:drizzlepac.astrodrizzle:	wcskey :	
INFO:drizzlepac.util:==== Processing Step  Static Mask  started at  01:12:53.708 (18/12/2019)
==== Processing Step  Static Mask  started at  01:12:53.708 (18/12/2019)
INFO:drizzlepac.staticMask:USER INPUT PARAMETERS for Static Mask Step:
INFO:drizzlepac.staticMask:	static :	True
INFO:drizzlepac.staticMask:	static_sig :	4.0
INFO:drizzlepac.staticMask:Computing static mask:

INFO:drizzlepac.staticMask:  mode = 55.075737;   rms = 10.970817;   static_sig = 4.00
INFO:drizzlepac.staticMask:  mode = 54.523428;   rms = 9.846490;   static_sig = 4.00
INFO:drizzlepac.staticMask:Computing static mask:

INFO:drizzlepac.staticMask:  mode = 50.450047;   rms = 10.694613;   static_sig = 4.00
INFO:drizzlepac.staticMask:  mode = 49.525925;   rms = 9.588069;   static_sig = 4.00
INFO:drizzlepac.staticMask:Computing static mask:

INFO:drizzlepac.staticMask:  mode = 51.066159;   rms = 10.823057;   static_sig = 4.00
INFO:drizzlepac.staticMask:  mode = 50.308009;   rms = 9.621310;   static_sig = 4.00
INFO:drizzlepac.staticMask:Computing static mask:

INFO:drizzlepac.staticMask:  mode = 50.493014;   rms = 11.489667;   static_sig = 4.00
INFO:drizzlepac.staticMask:  mode = 49.916463;   rms = 9.758092;   static_sig = 4.00
INFO:drizzlepac.staticMask:Saving static mask to disk: ./ACSWFC_2048x4096_1_staticMask.fits
INFO:drizzlepac.staticMask:Saving static mask to disk: ./ACSWFC_2048x4096_2_staticMask.fits
INFO:drizzlepac.util:==== Processing Step  Static Mask  finished at  01:12:55.31 (18/12/2019)
==== Processing Step  Static Mask  finished at  01:12:55.31 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Subtract Sky  started at  01:12:55.321 (18/12/2019)
==== Processing Step  Subtract Sky  started at  01:12:55.321 (18/12/2019)
INFO:drizzlepac.sky:USER INPUT PARAMETERS for Sky Subtraction Step:
INFO:drizzlepac.sky:	sky_bits :	0
INFO:drizzlepac.sky:	skyclip :	5
INFO:drizzlepac.sky:	skyfile :	
INFO:drizzlepac.sky:	skylower :	None
INFO:drizzlepac.sky:	skylsigma :	4.0
INFO:drizzlepac.sky:	skymask_cat :	
INFO:drizzlepac.sky:	skymethod :	localmin
INFO:drizzlepac.sky:	skystat :	median
INFO:drizzlepac.sky:	skysub :	True
INFO:drizzlepac.sky:	skyupper :	None
INFO:drizzlepac.sky:	skyuser :	
INFO:drizzlepac.sky:	skyusigma :	4.0
INFO:drizzlepac.sky:	skywidth :	0.1
INFO:drizzlepac.sky:	use_static :	True
INFO:stsci.skypac.utils:***** skymatch started on 2019-12-18 01:12:55.888421
***** skymatch started on 2019-12-18 01:12:55.888421
INFO:stsci.skypac.utils:      Version 1.0.4 (2019-05-30 14:31:23 -0400)
      Version 1.0.4 (2019-05-30 14:31:23 -0400)
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:'skymatch' task will apply computed sky differences to input image file(s).
'skymatch' task will apply computed sky differences to input image file(s).
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
INFO:stsci.skypac.utils:'MDRIZSKY' header keyword will represent sky value *computed* from data.
'MDRIZSKY' header keyword will represent sky value *computed* from data.
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:-----  User specified keywords:  -----
-----  User specified keywords:  -----
INFO:stsci.skypac.utils:       Sky Value Keyword:  'MDRIZSKY'
       Sky Value Keyword:  'MDRIZSKY'
INFO:stsci.skypac.utils:       Data Units Keyword: 'BUNIT'
       Data Units Keyword: 'BUNIT'
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:-----  Input file list:  -----
-----  Input file list:  -----
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   **  Input image: 'j8xi0xs0q_flc.fits'
   **  Input image: 'j8xi0xs0q_flc.fits'
INFO:stsci.skypac.utils:       EXT: 'SCI',1;	MASK: j8xi0xs0q_skymatch_mask_sci1.fits[0]
       EXT: 'SCI',1;	MASK: j8xi0xs0q_skymatch_mask_sci1.fits[0]
INFO:stsci.skypac.utils:       EXT: 'SCI',2;	MASK: j8xi0xs0q_skymatch_mask_sci2.fits[0]
       EXT: 'SCI',2;	MASK: j8xi0xs0q_skymatch_mask_sci2.fits[0]
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   **  Input image: 'j8xi0xs3q_flc.fits'
   **  Input image: 'j8xi0xs3q_flc.fits'
INFO:stsci.skypac.utils:       EXT: 'SCI',1;	MASK: j8xi0xs3q_skymatch_mask_sci1.fits[0]
       EXT: 'SCI',1;	MASK: j8xi0xs3q_skymatch_mask_sci1.fits[0]
INFO:stsci.skypac.utils:       EXT: 'SCI',2;	MASK: j8xi0xs3q_skymatch_mask_sci2.fits[0]
       EXT: 'SCI',2;	MASK: j8xi0xs3q_skymatch_mask_sci2.fits[0]
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   **  Input image: 'j8xi0xs6q_flc.fits'
   **  Input image: 'j8xi0xs6q_flc.fits'
INFO:stsci.skypac.utils:       EXT: 'SCI',1;	MASK: j8xi0xs6q_skymatch_mask_sci1.fits[0]
       EXT: 'SCI',1;	MASK: j8xi0xs6q_skymatch_mask_sci1.fits[0]
INFO:stsci.skypac.utils:       EXT: 'SCI',2;	MASK: j8xi0xs6q_skymatch_mask_sci2.fits[0]
       EXT: 'SCI',2;	MASK: j8xi0xs6q_skymatch_mask_sci2.fits[0]
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   **  Input image: 'j8xi0xsaq_flc.fits'
   **  Input image: 'j8xi0xsaq_flc.fits'
INFO:stsci.skypac.utils:       EXT: 'SCI',1;	MASK: j8xi0xsaq_skymatch_mask_sci1.fits[0]
       EXT: 'SCI',1;	MASK: j8xi0xsaq_skymatch_mask_sci1.fits[0]
INFO:stsci.skypac.utils:       EXT: 'SCI',2;	MASK: j8xi0xsaq_skymatch_mask_sci2.fits[0]
       EXT: 'SCI',2;	MASK: j8xi0xsaq_skymatch_mask_sci2.fits[0]
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:-----  Sky statistics parameters:  -----
-----  Sky statistics parameters:  -----
INFO:stsci.skypac.utils:       statistics function: 'median'
       statistics function: 'median'
INFO:stsci.skypac.utils:       lower = None
       lower = None
INFO:stsci.skypac.utils:       upper = None
       upper = None
INFO:stsci.skypac.utils:       nclip = 5
       nclip = 5
INFO:stsci.skypac.utils:       lsigma = 4.0
       lsigma = 4.0
INFO:stsci.skypac.utils:       usigma = 4.0
       usigma = 4.0
INFO:stsci.skypac.utils:       binwidth = 0.1
       binwidth = 0.1
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:-----  Data->Brightness conversion parameters for input files:  -----
-----  Data->Brightness conversion parameters for input files:  -----
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   *   Image: j8xi0xs0q_flc.fits
   *   Image: j8xi0xs0q_flc.fits
INFO:stsci.skypac.utils:       EXT = 'SCI',1
       EXT = 'SCI',1
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:       EXT = 'SCI',2
       EXT = 'SCI',2
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   *   Image: j8xi0xs3q_flc.fits
   *   Image: j8xi0xs3q_flc.fits
INFO:stsci.skypac.utils:       EXT = 'SCI',1
       EXT = 'SCI',1
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:       EXT = 'SCI',2
       EXT = 'SCI',2
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   *   Image: j8xi0xs6q_flc.fits
   *   Image: j8xi0xs6q_flc.fits
INFO:stsci.skypac.utils:       EXT = 'SCI',1
       EXT = 'SCI',1
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:       EXT = 'SCI',2
       EXT = 'SCI',2
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   *   Image: j8xi0xsaq_flc.fits
   *   Image: j8xi0xsaq_flc.fits
INFO:stsci.skypac.utils:       EXT = 'SCI',1
       EXT = 'SCI',1
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:       EXT = 'SCI',2
       EXT = 'SCI',2
INFO:stsci.skypac.utils:             Data units type: COUNTS
             Data units type: COUNTS
INFO:stsci.skypac.utils:             EXPTIME: 507.0 [s]
             EXPTIME: 507.0 [s]
INFO:stsci.skypac.utils:             Conversion factor (data->brightness):  0.7889546351084812
             Conversion factor (data->brightness):  0.7889546351084812
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:-----  Computing sky values requested image extensions (detector chips):  -----
-----  Computing sky values requested image extensions (detector chips):  -----
INFO:stsci.skypac.utils:

INFO:stsci.skypac.utils:   *   Image:   'j8xi0xs0q_flc.fits['SCI',1,2]'  --  SKY = 43.80914511294995 (brightness units)
   *   Image:   'j8xi0xs0q_flc.fits['SCI',1,2]'  --  SKY = 43.80914511294995 (brightness units)
INFO:stsci.skypac.utils:       Sky change (data units):
       Sky change (data units):
INFO:stsci.skypac.utils:      - EXT = 'SCI',1   delta(MDRIZSKY) = 55.5281   NEW MDRIZSKY = 55.5281
      - EXT = 'SCI',1   delta(MDRIZSKY) = 55.5281   NEW MDRIZSKY = 55.5281
INFO:stsci.skypac.utils:      - EXT = 'SCI',2   delta(MDRIZSKY) = 55.5281   NEW MDRIZSKY = 55.5281
      - EXT = 'SCI',2   delta(MDRIZSKY) = 55.5281   NEW MDRIZSKY = 55.5281
INFO:stsci.skypac.utils:   *   Image:   'j8xi0xs3q_flc.fits['SCI',1,2]'  --  SKY = 40.00184249125523 (brightness units)
   *   Image:   'j8xi0xs3q_flc.fits['SCI',1,2]'  --  SKY = 40.00184249125523 (brightness units)
INFO:stsci.skypac.utils:       Sky change (data units):
       Sky change (data units):
INFO:stsci.skypac.utils:      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.7023   NEW MDRIZSKY = 50.7023
      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.7023   NEW MDRIZSKY = 50.7023
INFO:stsci.skypac.utils:      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.7023   NEW MDRIZSKY = 50.7023
      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.7023   NEW MDRIZSKY = 50.7023
INFO:stsci.skypac.utils:   *   Image:   'j8xi0xs6q_flc.fits['SCI',1,2]'  --  SKY = 40.1377811469505 (brightness units)
   *   Image:   'j8xi0xs6q_flc.fits['SCI',1,2]'  --  SKY = 40.1377811469505 (brightness units)
INFO:stsci.skypac.utils:       Sky change (data units):
       Sky change (data units):
INFO:stsci.skypac.utils:      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.8746   NEW MDRIZSKY = 50.8746
      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.8746   NEW MDRIZSKY = 50.8746
INFO:stsci.skypac.utils:      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.8746   NEW MDRIZSKY = 50.8746
      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.8746   NEW MDRIZSKY = 50.8746
INFO:stsci.skypac.utils:   *   Image:   'j8xi0xsaq_flc.fits['SCI',1,2]'  --  SKY = 39.6654027453541 (brightness units)
   *   Image:   'j8xi0xsaq_flc.fits['SCI',1,2]'  --  SKY = 39.6654027453541 (brightness units)
INFO:stsci.skypac.utils:       Sky change (data units):
       Sky change (data units):
INFO:stsci.skypac.utils:      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.2759   NEW MDRIZSKY = 50.2759
      - EXT = 'SCI',1   delta(MDRIZSKY) = 50.2759   NEW MDRIZSKY = 50.2759
INFO:stsci.skypac.utils:      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.2759   NEW MDRIZSKY = 50.2759
      - EXT = 'SCI',2   delta(MDRIZSKY) = 50.2759   NEW MDRIZSKY = 50.2759
INFO:stsci.skypac.utils:***** skymatch ended on 2019-12-18 01:12:58.716886
***** skymatch ended on 2019-12-18 01:12:58.716886
INFO:stsci.skypac.utils:TOTAL RUN TIME: 0:00:02.828465
TOTAL RUN TIME: 0:00:02.828465
INFO:drizzlepac.util:==== Processing Step  Subtract Sky  finished at  01:12:59.345 (18/12/2019)
==== Processing Step  Subtract Sky  finished at  01:12:59.345 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Separate Drizzle  started at  01:12:59.347 (18/12/2019)
==== Processing Step  Separate Drizzle  started at  01:12:59.347 (18/12/2019)
INFO:drizzlepac.adrizzle:Interpreted paramDict with single=True as:
{'build': False, 'stepsize': 10, 'coeffs': True, 'wcskey': '', 'kernel': 'turbo', 'wt_scl': 'exptime', 'pixfrac': 1.0, 'fillval': None, 'bits': 336, 'compress': False, 'units': 'cps'}
INFO:drizzlepac.adrizzle:USER INPUT PARAMETERS for Separate Drizzle Step:
INFO:drizzlepac.adrizzle:	bits :	336
INFO:drizzlepac.adrizzle:	build :	False
INFO:drizzlepac.adrizzle:	clean :	True
INFO:drizzlepac.adrizzle:	coeffs :	True
INFO:drizzlepac.adrizzle:	compress :	False
INFO:drizzlepac.adrizzle:	crbit :	None
INFO:drizzlepac.adrizzle:	fillval :	None
INFO:drizzlepac.adrizzle:	kernel :	turbo
INFO:drizzlepac.adrizzle:	num_cores :	None
INFO:drizzlepac.adrizzle:	pixfrac :	1.0
INFO:drizzlepac.adrizzle:	proc_unit :	electrons
INFO:drizzlepac.adrizzle:	stepsize :	10
INFO:drizzlepac.adrizzle:	units :	cps
INFO:drizzlepac.adrizzle:	wcskey :	
INFO:drizzlepac.adrizzle:	wht_type :	None
INFO:drizzlepac.adrizzle:	wt_scl :	exptime
INFO:drizzlepac.adrizzle:  **Using sub-sampling value of 10 for kernel turbo
INFO:drizzlepac.adrizzle:Running Drizzle to create output frame with WCS of: 
INFO:astropy.wcs.wcs:WCS Keywords
WCS Keywords
INFO:astropy.wcs.wcs:

INFO:astropy.wcs.wcs:Number of WCS axes: 2
Number of WCS axes: 2
INFO:astropy.wcs.wcs:CTYPE : 'RA---TAN'  'DEC--TAN'  
CTYPE : 'RA---TAN'  'DEC--TAN'  
INFO:astropy.wcs.wcs:CRVAL : 150.20204715246223  2.6769864498720612  
CRVAL : 150.20204715246223  2.6769864498720612  
INFO:astropy.wcs.wcs:CRPIX : 2111.5  2210.0  
CRPIX : 2111.5  2210.0  
INFO:astropy.wcs.wcs:CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
INFO:astropy.wcs.wcs:CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
INFO:astropy.wcs.wcs:NAXIS : 4223  4420
NAXIS : 4223  4420
INFO:drizzlepac.adrizzle:Executing 4 parallel workers
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 55.528091 to j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.702335 to j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.874638 to j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.275898 to j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.275898 to j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.874638 to j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 55.528091 to j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.702335 to j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.cdriz:-Drizzling using kernel = turbo

INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs6q_single_sci.fits
-Generating simple FITS output: j8xi0xs6q_single_sci.fits
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs0q_single_sci.fits
-Generating simple FITS output: j8xi0xs0q_single_sci.fits
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xsaq_single_sci.fits
-Generating simple FITS output: j8xi0xsaq_single_sci.fits
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs3q_single_sci.fits
-Generating simple FITS output: j8xi0xs3q_single_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs6q_single_sci.fits
Writing out image to disk: j8xi0xs6q_single_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs0q_single_sci.fits
Writing out image to disk: j8xi0xs0q_single_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xsaq_single_sci.fits
Writing out image to disk: j8xi0xsaq_single_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs3q_single_sci.fits
Writing out image to disk: j8xi0xs3q_single_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs6q_single_wht.fits
Writing out image to disk: j8xi0xs6q_single_wht.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs0q_single_wht.fits
Writing out image to disk: j8xi0xs0q_single_wht.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xsaq_single_wht.fits
Writing out image to disk: j8xi0xsaq_single_wht.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs3q_single_wht.fits
Writing out image to disk: j8xi0xs3q_single_wht.fits
INFO:drizzlepac.util:==== Processing Step  Separate Drizzle  finished at  01:13:03.214 (18/12/2019)
==== Processing Step  Separate Drizzle  finished at  01:13:03.214 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Create Median  started at  01:13:03.220 (18/12/2019)
==== Processing Step  Create Median  started at  01:13:03.220 (18/12/2019)
INFO:drizzlepac.createMedian:USER INPUT PARAMETERS for Create Median Step:
INFO:drizzlepac.createMedian:	combine_bufsize :	None
INFO:drizzlepac.createMedian:	combine_grow :	1
INFO:drizzlepac.createMedian:	combine_hthresh :	None
INFO:drizzlepac.createMedian:	combine_lthresh :	None
INFO:drizzlepac.createMedian:	combine_maskpt :	0.3
INFO:drizzlepac.createMedian:	combine_nhigh :	0
INFO:drizzlepac.createMedian:	combine_nlow :	0
INFO:drizzlepac.createMedian:	combine_nsigma :	4 3
INFO:drizzlepac.createMedian:	combine_type :	median
INFO:drizzlepac.createMedian:	compress :	False
INFO:drizzlepac.createMedian:	median :	True
INFO:drizzlepac.createMedian:	median_newmasks :	True
INFO:drizzlepac.createMedian:	proc_unit :	native
INFO:drizzlepac.createMedian:reference sky value for image 'j8xi0xs0q_flc.fits' is 55.52809143066406
reference sky value for image 'j8xi0xs0q_flc.fits' is 55.52809143066406
INFO:drizzlepac.createMedian:reference sky value for image 'j8xi0xs3q_flc.fits' is 50.70233535766601
reference sky value for image 'j8xi0xs3q_flc.fits' is 50.70233535766601
INFO:drizzlepac.createMedian:reference sky value for image 'j8xi0xs6q_flc.fits' is 50.874637603759766
reference sky value for image 'j8xi0xs6q_flc.fits' is 50.874637603759766
INFO:drizzlepac.createMedian:reference sky value for image 'j8xi0xsaq_flc.fits' is 50.27589797973633
reference sky value for image 'j8xi0xsaq_flc.fits' is 50.27589797973633
INFO:drizzlepac.createMedian:Saving output median image to: 'f814w_med.fits'
Saving output median image to: 'f814w_med.fits'
INFO:drizzlepac.util:==== Processing Step  Create Median  finished at  01:13:10.93 (18/12/2019)
==== Processing Step  Create Median  finished at  01:13:10.93 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Blot  started at  01:13:10.933 (18/12/2019)
==== Processing Step  Blot  started at  01:13:10.933 (18/12/2019)
INFO:drizzlepac.ablot:USER INPUT PARAMETERS for Blot Step:
INFO:drizzlepac.ablot:	blot_addsky :	True
INFO:drizzlepac.ablot:	blot_interp :	poly5
INFO:drizzlepac.ablot:	blot_sinscl :	1.0
INFO:drizzlepac.ablot:	blot_skyval :	0.0
INFO:drizzlepac.ablot:	coeffs :	True
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs0q_flc.fits[sci,1]
    Blot: creating blotted image:  j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 55.528091 to blotted image j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs0q_sci1_blt.fits
-Generating simple FITS output: j8xi0xs0q_sci1_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs0q_sci1_blt.fits
Writing out image to disk: j8xi0xs0q_sci1_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs0q_flc.fits[sci,2]
    Blot: creating blotted image:  j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 55.528091 to blotted image j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs0q_sci2_blt.fits
-Generating simple FITS output: j8xi0xs0q_sci2_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs0q_sci2_blt.fits
Writing out image to disk: j8xi0xs0q_sci2_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs3q_flc.fits[sci,1]
    Blot: creating blotted image:  j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.702335 to blotted image j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs3q_sci1_blt.fits
-Generating simple FITS output: j8xi0xs3q_sci1_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs3q_sci1_blt.fits
Writing out image to disk: j8xi0xs3q_sci1_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs3q_flc.fits[sci,2]
    Blot: creating blotted image:  j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.702335 to blotted image j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs3q_sci2_blt.fits
-Generating simple FITS output: j8xi0xs3q_sci2_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs3q_sci2_blt.fits
Writing out image to disk: j8xi0xs3q_sci2_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs6q_flc.fits[sci,1]
    Blot: creating blotted image:  j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.874638 to blotted image j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs6q_sci1_blt.fits
-Generating simple FITS output: j8xi0xs6q_sci1_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs6q_sci1_blt.fits
Writing out image to disk: j8xi0xs6q_sci1_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xs6q_flc.fits[sci,2]
    Blot: creating blotted image:  j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.874638 to blotted image j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xs6q_sci2_blt.fits
-Generating simple FITS output: j8xi0xs6q_sci2_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xs6q_sci2_blt.fits
Writing out image to disk: j8xi0xs6q_sci2_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xsaq_flc.fits[sci,1]
    Blot: creating blotted image:  j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.275898 to blotted image j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xsaq_sci1_blt.fits
-Generating simple FITS output: j8xi0xsaq_sci1_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xsaq_sci1_blt.fits
Writing out image to disk: j8xi0xsaq_sci1_blt.fits
INFO:drizzlepac.ablot:    Blot: creating blotted image:  j8xi0xsaq_flc.fits[sci,2]
    Blot: creating blotted image:  j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.ablot:Using default C-based coordinate transformation...
Using default C-based coordinate transformation...
INFO:drizzlepac.ablot:Applying sky value of 50.275898 to blotted image j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.outputimage:-Generating simple FITS output: j8xi0xsaq_sci2_blt.fits
-Generating simple FITS output: j8xi0xsaq_sci2_blt.fits
INFO:drizzlepac.outputimage:Writing out image to disk: j8xi0xsaq_sci2_blt.fits
Writing out image to disk: j8xi0xsaq_sci2_blt.fits
INFO:drizzlepac.util:==== Processing Step  Blot  finished at  01:13:22.826 (18/12/2019)
==== Processing Step  Blot  finished at  01:13:22.826 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Driz_CR  started at  01:13:22.829 (18/12/2019)
==== Processing Step  Driz_CR  started at  01:13:22.829 (18/12/2019)
INFO:drizzlepac.drizCR:USER INPUT PARAMETERS for Driz_CR Step:
INFO:drizzlepac.drizCR:	crbit :	4096
INFO:drizzlepac.drizCR:	driz_cr :	True
INFO:drizzlepac.drizCR:	driz_cr_corr :	False
INFO:drizzlepac.drizCR:	driz_cr_ctegrow :	0
INFO:drizzlepac.drizCR:	driz_cr_grow :	1
INFO:drizzlepac.drizCR:	driz_cr_scale :	1.2 0.7
INFO:drizzlepac.drizCR:	driz_cr_snr :	3.5 3.0
INFO:drizzlepac.drizCR:	inmemory :	False
INFO:drizzlepac.drizCR:Executing 4 parallel workers
INFO:drizzlepac.drizCR:Creating output: j8xi0xs0q_sci1_crmask.fits
Creating output: j8xi0xs0q_sci1_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xsaq_sci1_crmask.fits
Creating output: j8xi0xsaq_sci1_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xs6q_sci1_crmask.fits
Creating output: j8xi0xs6q_sci1_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xs3q_sci1_crmask.fits
Creating output: j8xi0xs3q_sci1_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xsaq_sci2_crmask.fits
Creating output: j8xi0xsaq_sci2_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xs0q_sci2_crmask.fits
Creating output: j8xi0xs0q_sci2_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xs3q_sci2_crmask.fits
Creating output: j8xi0xs3q_sci2_crmask.fits
INFO:drizzlepac.drizCR:Creating output: j8xi0xs6q_sci2_crmask.fits
Creating output: j8xi0xs6q_sci2_crmask.fits
INFO:drizzlepac.util:==== Processing Step  Driz_CR  finished at  01:13:27.982 (18/12/2019)
==== Processing Step  Driz_CR  finished at  01:13:27.982 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.util:==== Processing Step  Final Drizzle  started at  01:13:27.993 (18/12/2019)
==== Processing Step  Final Drizzle  started at  01:13:27.993 (18/12/2019)
INFO:drizzlepac.adrizzle:Interpreted paramDict with single=False as:
{'build': False, 'stepsize': 10, 'coeffs': True, 'wcskey': '', 'wht_type': 'EXP', 'kernel': 'square', 'wt_scl': 'exptime', 'pixfrac': 1.0, 'fillval': None, 'maskval': None, 'bits': 336, 'units': 'cps'}
INFO:drizzlepac.adrizzle:USER INPUT PARAMETERS for Final Drizzle Step:
INFO:drizzlepac.adrizzle:	bits :	336
INFO:drizzlepac.adrizzle:	build :	False
INFO:drizzlepac.adrizzle:	clean :	True
INFO:drizzlepac.adrizzle:	coeffs :	True
INFO:drizzlepac.adrizzle:	crbit :	4096
INFO:drizzlepac.adrizzle:	fillval :	None
INFO:drizzlepac.adrizzle:	kernel :	square
INFO:drizzlepac.adrizzle:	maskval :	None
INFO:drizzlepac.adrizzle:	pixfrac :	1.0
INFO:drizzlepac.adrizzle:	proc_unit :	native
INFO:drizzlepac.adrizzle:	stepsize :	10
INFO:drizzlepac.adrizzle:	units :	cps
INFO:drizzlepac.adrizzle:	wcskey :	
INFO:drizzlepac.adrizzle:	wht_type :	EXP
INFO:drizzlepac.adrizzle:	wt_scl :	exptime
INFO:drizzlepac.adrizzle:  **Using sub-sampling value of 10 for kernel square
INFO:drizzlepac.adrizzle:Running Drizzle to create output frame with WCS of: 
INFO:astropy.wcs.wcs:WCS Keywords
WCS Keywords
INFO:astropy.wcs.wcs:

INFO:astropy.wcs.wcs:Number of WCS axes: 2
Number of WCS axes: 2
INFO:astropy.wcs.wcs:CTYPE : 'RA---TAN'  'DEC--TAN'  
CTYPE : 'RA---TAN'  'DEC--TAN'  
INFO:astropy.wcs.wcs:CRVAL : 150.20204715246223  2.6769864498720612  
CRVAL : 150.20204715246223  2.6769864498720612  
INFO:astropy.wcs.wcs:CRPIX : 2111.5  2210.0  
CRPIX : 2111.5  2210.0  
INFO:astropy.wcs.wcs:CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
CD1_1 CD1_2  : 2.4116113708575043e-06  1.3677915234561588e-05  
INFO:astropy.wcs.wcs:CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
CD2_1 CD2_2  : 1.3677915234561588e-05  -2.4116113708575043e-06  
INFO:astropy.wcs.wcs:NAXIS : 4223  4420
NAXIS : 4223  4420
INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 55.528091 to j8xi0xs0q_flc.fits[sci,1]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 1
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 55.528091 to j8xi0xs0q_flc.fits[sci,2]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 2
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.702335 to j8xi0xs3q_flc.fits[sci,1]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 1
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.702335 to j8xi0xs3q_flc.fits[sci,2]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 2
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.874638 to j8xi0xs6q_flc.fits[sci,1]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 1
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.874638 to j8xi0xs6q_flc.fits[sci,2]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 2
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.adrizzle:Applying sky value of 50.275898 to j8xi0xsaq_flc.fits[sci,1]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 1
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.adrizzle:-Drizzle input: j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.adrizzle:Applying sky value of 50.275898 to j8xi0xsaq_flc.fits[sci,2]
INFO:drizzlepac.imageObject:Applying EXPTIME weighting to DQ mask for chip 2
INFO:drizzlepac.adrizzle:Using WCSLIB-based coordinate transformation...
INFO:drizzlepac.adrizzle:stepsize = 10
INFO:drizzlepac.cdriz:-Drizzling using kernel = square

INFO:drizzlepac.outputimage:-Generating simple FITS output: f814w_drc_sci.fits
-Generating simple FITS output: f814w_drc_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: f814w_drc_sci.fits
Writing out image to disk: f814w_drc_sci.fits
INFO:drizzlepac.outputimage:Writing out image to disk: f814w_drc_wht.fits
Writing out image to disk: f814w_drc_wht.fits
INFO:drizzlepac.util:==== Processing Step  Final Drizzle  finished at  01:13:50.689 (18/12/2019)
==== Processing Step  Final Drizzle  finished at  01:13:50.689 (18/12/2019)
INFO:drizzlepac.util:

INFO:drizzlepac.astrodrizzle:

INFO:drizzlepac.astrodrizzle:AstroDrizzle Version 3.1.3 is finished processing at 01:13:50.693 (18/12/2019).
AstroDrizzle Version 3.1.3 is finished processing at 01:13:50.693 (18/12/2019).
INFO:drizzlepac.astrodrizzle:

INFO:drizzlepac.util:

INFO:drizzlepac.util:   --------------------          --------------------
   --------------------          --------------------
INFO:drizzlepac.util:                   Step          Elapsed time
                   Step          Elapsed time
INFO:drizzlepac.util:   --------------------          --------------------
   --------------------          --------------------
INFO:drizzlepac.util:

INFO:drizzlepac.util:         Initialization          1.8117 sec.
         Initialization          1.8117 sec.
INFO:drizzlepac.util:            Static Mask          1.6103 sec.
            Static Mask          1.6103 sec.
INFO:drizzlepac.util:           Subtract Sky          4.0237 sec.
           Subtract Sky          4.0237 sec.
INFO:drizzlepac.util:       Separate Drizzle          3.8669 sec.
       Separate Drizzle          3.8669 sec.
INFO:drizzlepac.util:          Create Median          7.7116 sec.
          Create Median          7.7116 sec.
INFO:drizzlepac.util:                   Blot          11.8928 sec.
                   Blot          11.8928 sec.
INFO:drizzlepac.util:                Driz_CR          5.1534 sec.
                Driz_CR          5.1534 sec.
INFO:drizzlepac.util:          Final Drizzle          22.6961 sec.
          Final Drizzle          22.6961 sec.
INFO:drizzlepac.util:   ====================          ====================
   ====================          ====================
INFO:drizzlepac.util:                  Total          58.7665 sec.
                  Total          58.7665 sec.
INFO:drizzlepac.imageObject:Removing intermediate files for j8xi0xs0q_flc.fits
INFO:drizzlepac.imageObject:Removing intermediate files for j8xi0xs3q_flc.fits
INFO:drizzlepac.imageObject:Removing intermediate files for j8xi0xs6q_flc.fits
INFO:drizzlepac.imageObject:Removing intermediate files for j8xi0xsaq_flc.fits
INFO:drizzlepac.util:Trailer file written to:  f814w_driz.log
Trailer file written to:  f814w_driz.log

5. Display the original FLC image and the combined DRC data

In [24]:
sci2 = fits.getdata('j8xi0xsaq_flc.fits', ('sci', 2))
sci1 = fits.getdata('j8xi0xsaq_flc.fits', ('sci', 1))

fig = plt.figure(figsize=(50, 50))
ax1 = fig.add_subplot(2, 2, 1)
ax2 = fig.add_subplot(1, 2, 1)

ax1.imshow(sci2, vmin=50, vmax=100, cmap='Greys_r', origin='lower')
ax2.imshow(sci1, vmin=50, vmax=100, cmap='Greys_r', origin='lower')
Out[24]:
<matplotlib.image.AxesImage at 0x7f8be1b0c510>
In [25]:
sci = fits.getdata('f814w_drc_sci.fits')
fig = plt.figure(figsize=(20, 20))
plt.imshow(sci, vmin=-0.01, vmax=0.1, cmap='Greys_r', origin='lower')
Out[25]:
<matplotlib.image.AxesImage at 0x7f8c0328e6d0>

Summary

TweakReg may be used to align HST images based on the position of objects in the frame. Point sources are typically used for this purpose, but compact sources such as background galaxies may also be used by increasing the value of the parameter conv_width in imagefindpars. The data quality arrays of the input calibrated frames may also be used to further improve the fits by telling TweakReg to ignore pixels with specific flags inthe DQ array via the parameter dqbits in imagefindpars.

About this Notebook

Author: J. Mack, STScI WFC3 Team
Updated: January 8, 2019